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,— , Abstract 



We present new classical generalized one-monopole solution of the SU(2) 
Yang-Mills-Higgs theory with the Higgs field in the adjoint representation. 
We show that this generalized solution with 0- winding number m = 1 and (p- 
winding number n = 1 is an axially symmetric Jacobi elliptic generalization 
, of the 't Hooft-Polyakov one-monopole. We construct this axially symmetric 

^ one-monopole solution by generalizing the large distance asymptotic solution 

of the 't Hooft-Polyakov one-monopole to the Jacobi elliptic functions and 
solving the second order equations of motion numerically when the Higgs 
■nj- potential is vanishing and non vanishing. These solutions are regular non- 

BPS finite energy solutions. 

o 
o 

Tl 1 Introduction 

> 

^ The 't Hooft-Polyakov one-monopole solution of the SU(2) Yang-Mills-Higgs (YMH) 

field theory with nonzero Higgs mass and self-interaction is the first monopole so- 
lution to possess finite energy [T]. This numerical, spherically symmetric non-BPS 
one monopole solution is invariant under a U(l) subgroup of the local SU(2) gauge 
group. All known finite energy solutions of the YMH field theory with a unit mag- 
netic charge are spherically symmetric [l]-[2]. All other finite energy monopole 
configurations with magnetic charges greater than unity do not possess spherical 
symmetry [3] as it has been shown that these solutions cannot possess spherical 
symmetry when energy is finite 

Exact solutions of the YMH equations only exist in the Bogomol'nyi-Prasad- 
Sommerfield (BPS) limit when the Higgs potential vanishes or when the Higgs 
field self-interaction constant A = 0. Outside the BPS limit and when A 7^ 0, only 
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numerical solutions exist. The monopole-antimonopole pair (MAP), monopole- 
antimonopole chain (MAC) and vortex rings configurations of Ref. [5] are non- 
BogomoPnyi solutions and they are numerical solutions even in the limit of van- 
ishing Higgs potential. 

Here, we present new classical generalized Jacobi elliptic axially symmetric 
one-monopole solutions. These new monopoles are regular solutions of the SU(2) 
YMH theory with finite energies. Similar to the monopole solutions of Ref. [5], 
they are non-Bogomol'nyi solutions and are numerical solutions even in the limit 
of vanishing Higgs potential. 

2 The SU(2) Yang- Mills- Higgs Theory 

The SU(2) YMH Lagrangian in 3+1 dimensions with non vanishing Higgs potential 
is given by 



£ = _ - iA(<i>"<i>" - ^)'\ (1) 

4 2 4 A 



Here the Higgs field mass is fi and the strength of the Higgs potential is A which 
are constants. The vacuum expectation value of the Higgs field is ^ = /i/a/A. 
The Lagrangian ([T]) is gauge invariant under the set of independent local SU(2) 
transformations at each space-time point. The covariant derivative of the Higgs 
field and the gauge field strength tensor are given respectively by 

F;, = d,A:-d.A; + ge^'^AlAl. (2) 

Since the gauge field coupling constant g can be scaled away, we can set g to 
one without any loss of generality. The metric used is = ( — h ++). The 
SU(2) internal group indices a,b,c run from 1 to 3 and the spatial indices are 
^.v.a = 0, 1, 2, and 3 in Minkowski space. 

The equations of motion that follow from the Lagrangian ([T]) are 

,,2 

= A$"(<l>''<l>^- ^). (3) 

A 

Non-BPS solutions to the YMH theory are obtained by solving the second order 
differential equations of motion ([3]), whereas BPS solutions can be obtained more 
easily by solving the first order Bogomol'nyi equations, 

Bf ± Di^" = 0. (4) 

The static energy for the system is 



2\ 2 



(5) 
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and in the BPS limit when the Higgs potential vanishes, 

£ = ^ j 9,(5^$'^) <fx + J ^{B^ ± Di<^>''f £x 

= tJ di{Bf^^) d^x = AttM^, (6) 

where A4 is the "topological charge" . The tensor identified with the electromag- 
netic field, as proposed by 't Hooft [2] is 

= d,A,-d,A^-e''''^''d^^%^', (7) 

where = ^''A^, the Higgs unit vector, = $"/|$|, and the Higgs field mag- 
nitude |$| = y/ The Abelian electric field is Ei = Fq^, and the Abelian 
magnetic field is Bi = —^eij^Fj^. We would also like to write the Abelian 't Hooft 
electromagnetic field as 

F^u = G^u + Hfj,^, (8) 

where 

G^u = d^Ajj — dyA^, 

H^, = -e'''"'¥d^^''d,^\ (9) 

which we refer to as the gauge part and the Higgs part of the 't Hooft electomag- 
netic field respectively. 

The topological magnetic current [6] is defined to be 

K = ^ ^abc d''^ d^^' 9^$^ (10) 

which is also the topological current density of the system and the corresponding 
conserved topological magnetic charge is 



= ^<fd'a,{e.,,,e'^''^^''d,^%^ 



''x 



As mentioned by Arafune et al. [7|, the magnetic charge Ai is the total magnetic 
charge of the system if and only if the gauge field is nonsingular. If the gauge field 
is singular and carries Dirac string monopoles, then the total magnetic charge of 
the system is the sum of the Dirac string monopoles and the monopoles carry by 
the Higgs field. 
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3 The Exact Generalized Asymptotic Solutions 



The magnetic ansatz of Ref. [5] possesses internal space unit vectors with ^-winding 
number m and 0-winding number n greater than one, 

Al = 0, = $i(r,0) {i^ + $2(r,^){i^, (12) 

where the spatial spherical coordinate orthonormal unit vectors are 

fj = sin 6' COS0 5ii + sin^ sin0 5i2 + cos^ (5^3, 
9i = cos 9 cos (p 6ii + cos 9 sin 5^2 — sin 6* 5^3, 
0i = - sin0 5ji + COS0 (5i2, (13) 

and the isospin coordinate orthonormal unit vectors are 

■u^ = sin m6' cos n(j) 61 + sin m9 sin n(j) 62 + cos 7719 5^, 
Ug = cos m9 cos ncj) 51 + cos m9 sinnc/) ^2 sinm^^ ^3, 

= — sinn0 51 + cosncj) where m > 1, n > 1. (14) 

In Ref.[S], the MAP solutions and the MAC solutions corresponds to MA- MA- - 

MA pairs and M-A-M-A- -M-A-M chain respectively lying on the 2;-axis. Hence 

the MAP solutions possess zero net magnetic charge whereas the MAC solutions 
have a net magnetic charge of n = 1, 2. When m = 2, 4, 6, the solutions are the 
1-MAP, 2-MAP, 3-MAP, configurations respectively and when m = 1, 3, 5, 

the solutions are the M, M-A-M, M- A- M-A-M, , monopole chains respectively. 

Kleihaus et al. [5] have shown that these MAP and MAC solutions can have a 
single pole charge of n = 1,2. When the 0-winding number n = 3 and above, 
vortex rings are formed when m > 2. 
By using the relationship [8j 

= cos(m - 1)9 hi + sin(m - 1)6* hi, 
itg = - sin(m - 1)9 h^. + cos(m - 1)^ n^, ul = hi (15) 



where n^, hg, and h'^ are the unit vectors of Eq.( 14) when m = 1, we can reduce the 
MAP and MAC solutions of Ref. to the 6'-winding number, m = 1, solutions. 
Hence the asymptotic solutions with m = 1 at large r is given by 



R2 



1 + ^2 



?2(sin^ + sin p6'(a cos ^ + h)) 



sin 6* 

n(cos 9 — cosp9{a cos 9 + b)) 
sin 9 ' 



, Ri = 0, 



= ^cosp9, ^2 = ^sinp9, (16) 



with a = 0, 6 = 1, and p = 1, 3, 5, ... in the MAP solution and a = 1, 6 = 0, and 
p = 0,2,4, ... in the MAC solution. Here we notice that the solution (16) can be 
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generalized to 

dgh{oo,6) n{sm6 + h{oo,6){acos6 + b)) 

¥^1 = 1 H 7 ^) W2 = ^ , 

g[oo,U) smu 

no ?^(cos — (7(00, ^^)(a COS 6* + 6)) 

Ri = 0, R2 = — , 

smu 

$1 = ^gioo,9), ^2=^h{oo,e), (17) 

where g{r, 6)^ + h{r, 9)'^ = 1 and a, b, ^ are constants. In Ref.jHj, we notice that the 
solution can be generahzed to include the Jacobi eUiptic functions Ea = Ea{u, k); 
A = 1,2, as follow, 

g{oo,e) = Ei{u,k) and h{oo,e) = E2{u,k). (18) 

Here u = p6 and Ei{u,k) and E2{u,k) are a pair of non-singular Jacobi elliptic 
functions that satisfy the relation Ei{u, k)"^ + E2{u, k)"^ = 1. 

The three possible pairs of non-singular Jacobi elliptic solutions are 

Ei{u,k) = cn{u,k) and E2{u, k) = sn{u, k) (19) 
Ei{u,k) = dn{u,k) and E2{u, k) = k sn{u, k) (20) 
^ , cn(u,k) , ^ , ,^ k'sniu.k) , 

where < < 1 is the Jacobi elliptic parameter and k' = a/1 — k^- We then label 
the three sets of solutions (19) - (21) by calling them the Jacobi elliptic A (JEA), 
Jacobi elliptic B (JEB), and Jacobi elliptic C (JEC) solutions respectively. One 
should also notice that when k = 0, the JEA and JEC solutions revert back to 
solution (16) and the JEB solution becomes the 't Hooft-Polyakov one-monopole 
solution when a = 1, 6 = 0, and n = 1. 

We have also notice that an even more general asymptotic solution exists when 

g{oo,9) = sin q9 E2{u,k) + cos q9 Ei{u,k), (22) 
h{oo,9) = cos q9 E2{u,k) — sin q9 Ei{u,k), (23) 



where g is a constant. When Ea{u, k) is given by Eq.(19), we have the JEA type 
solution. The gauge potentials and Higgs fields of the one-monopole JEA type 
solution at large r is then given by 

A" = --{l-q + pdn(p9,k)}4>''9i 
r 

+ -{sin6' + {cosg6' sn(p9,k) — sin q9 cnip9,k)\ cos 9} 9"' (pi 

r sm 9 

I {cos 6' — {sing6' sn{p9, k) + cos q9 cn{p9, k)} cos 9} f "'(pi, 



r sin 9 



^0 = 0, 



= ^{sin q9 sn{p9, k) + cos q9 cn{p9, k)} f" 

+ ^{cos q9 sn{p9,k) - sin q9 cn{p9,k)}9". (24) 
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By using suitable choices for the solution constants {p,q,k), the gauge potentials 
of Eq.(24) can be made non singular. In this paper we solved numerically for four 
one-monopole JEA type solutions when k = ki = 0.9844325133 or when the period 
of the Jacobi elliptic function A is An, 

/ ; = TT. 25 

^0 Vl - A;2 sin^ 6 

Four axially symmetric one-monopole solutions are obtained when (j9, g, k) = 
(1, |,fci),(l,§,fci),(2, l,fci), and (2,3,A;i). When {p,q,k) = (0,0, fc), the solu- 
tion is just the 't Hooft-Polyakov one-monopole. The solution (24) satisfied the 
temporal gauge condition, Aq = 0, only. The non-Abelian magnetic fields and the 
Higgs field covariant derivatives are given respectively by 

+ + (V-i - ^2) cot + (1 - ^i)R2)es-i 
+ + {cote- R2)Ri)e''e, + ^{rij[ + Ri)rk 

= Bi^"" (26) 
= {<!>[- ^Ri<^2Kn + {<l>'^ + ^Ri<^i)e^n 

+ ^($1 - (1 - ^i)<l>2)r% + h^2 + (1 - 

+ -{{1- ^2)^1 + {cote -R2)^2)r^i 

r 

= 0, (27) 
From Eq.(|7[ - (|9|), the 't Hooft Abelian magnetic field Bi is given by 
B, = Bf + Bf, 

Bf = -e,,kdjAk and B^ = ^e,,fce"'^<l'^9,<l''a,l>^ (28) 

In Ref.jS], it mentioned that is the Abelian component of A°^, and that it is 
the unrestricted part of the gauge potential A", that is, it is not restricted by the 
condition, 

= d^^" + e^'^'Al^" = 0. 



In fact the general large r asymptotic solution (17) can be written as 

where A^ = A (f)^, = g + h 

1 , , s. h sin e — q cos e + a cos e + h 

A = -{tfj, h - R2 g) = — . 

r r sm e 



The other component of that is restricted is — e'^^^<|)^c}^<|)^ and it is completely 
determined by the magnetic symmetry. Hence the unrestricted component of 
is referred to as the electric part and the restricted component of is referred to 
as the magnetic part of the gauge potential 0. 

With some calculations, it can be shown that the magnetic part of i?" which 
is also the Higgs part is 

Bf = —neijkdj {g cos 6 — hsin 9) 9^0, (29) 

and the electric part which is also the gauge part is 

Bf = -neijkdj{^^{i^2h-R2g)^ dkcj), (30) 

where Bf is from the magnetic component of and Bf is from the electric 
component of Therefore the net magnetic field when n = 1 is 

Bi = -neijkdj {cos a{r, 6)} dk4>, 
where cosQ;(r,6') = sin^(5f(cot 6* — i?2) — ^(1 — ^2))- (31) 

A contour plot of cos a{r, 9) = constant at a fixed value of will give the magnetic 
field lines of the one-monopole solutions. In the limit when r ^ 00, n = a = 1 
and 6 = 0, 

fi 



B^ = -eijkdj {cos 6} d, 



is just the magnetic field of a one-monopole. The Abelian energy density of the 
monopole is then given by 

^BiB' = ^ 2 n ^j cos a cos a. (32) 
2 2r2sm^^ 

4 The Numerical One-Monopole Solutions 

The numerical calculations were performed using the Maple 12 and MatLab R2009a 
softwares. The second order equations of motion dsj) are reduced to six partial dif- 



ferential equations with the ansatz (12) which are then transformed into a system 
of nonlinear equations using the finite difference approximation. This system of 
nonlinear equations are then discretized on a non-equidistant grid of size 70 x 60 
covering the integration regions < a; < 1 and < ^ < tt. Here x is the finite 
interval compactified coordinate given by a; = The partial derivative with re- 
spect to the radial coordinate is then replaced accordingly by dr — )■ (1 — xYd^ and 
^ ^ {1 — xY-§^ — 2(1 — xY-§^- First of all, we used Maple to find the Jacobian 
sparsity pattern for the system of nonlinear equations. After that we provide this 
information to Matlab to run the numerical computation. The system of nonlinear 
equations are then solved numerically using the trust-region-reflective algorithm 
by providing the solver with good initial guess. 
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z 




The second order equations of motion Eq. (|3j) were solved when the winding 
numbers m = 1 and n = 1 and firstly with zero Higgs potential, that is A = /i = 
but with nonzero expectation value ^ = 1 and finally with A = /i = 1 and ^ = 1. 
Hence the solutions obtained are non-BPS solutions with total energy given by 



The boundary conditions used at small r is given by 

^i(0,^) =V'2(0,e) = i?i(0,^) =i?2(0,e) 
$i(0,^) =^ocos^, $2(0,0) = -^0 sin ^ 
sine $i(O,0) + COS0 $2(0,0) 

d 



(33) 



dr 



{cos6 $i(r, 6) - sin^ ^2{r, 9)) 








(34) 



r=0 



and at large r is given by Eq.(24). In order for the solution to be regular along 
the z-axis, the constraint imposed on the solution is 



=0,7r 



=0,7r 



=0,7r 
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Air, 9) 



e=o,TT 



89 



^2(r, 9) 



d_ 
89 



(35) 



=0,7r 



The above constraints are sufficient to ensure the regularity of the solution and 
the gauge fixing condition 



r—Ri{r,9) - -^Mr,0) 



dr 



39 



0, 



(36) 



is not imposed upon the numerical solutions here as Eq.(36) will restrict the one- 



monopole solutions to just one solution, that is the 't Hooft-Polyakov one-monopole 



solution. Similarly if the gauge fixing condition Eq.(36) is imposed on the MAP 



and MAC solutions, then these solutions are restricted to those discussed in Ref. 
[5]. We also note that the MAP and MAC solutions of Ref. |5] does not satisfied 
the radiation gauge condition d^AI = in general. It is only the 't Hooft-Polyakov 
one-monopole and the one-MAP solution that will satisfied the radiation gauge 
when the winding number n = 1. All the other MAC and MAP solutions do not 
satisfied the radiation gauge condition. 



The boundaries conditions imposed on the gauge potentials at large r, Eq.(24) 



small r, Eq.(34), and along the 2;-axis, Eq.(35), are sufficient to ensure regularity of 



the solutions. The ansatz ( 12 ) will automatically ensure that the 't Hooft Abelian 
gauge potential. 



A- = 



1 



{'4,2{r,9)h{r,9) - R2{r,9)g{r,9))^, 



(37) 



obeys the radiation gauge condition, d'^Ai = 0. 

By using suitable initial values, we obtained the new one-monopole numerical 



solutions when the the constants (p, q, k) of the asymptotic solution (24) takes the 
following sets of values 



(p, q, k) 
ki 



^1) 2' 



^1) 2 ' 



(2,l,fci), (2,3,fci) 



0.9844325133. 



(38) 



These JEA one-monopole solutions are axially symmetric monopoles with Jacobi 
elliptic functions of period An in the asymptotic large r solutions. Upon solving 
these four solutions using our numerical procedures, we plotted the polar plot 
of the magnitude of the Higgs magnetic field \r'^Bf\ of these monopole solutions 
together with the 't Hooft-Polyakov monopole's \r^Bl^\ versus 9, at large r, Figjl} 
We also calculate for the energy density. 



ED{r, 9) = -{B^B^ + A$"A$") + -A($"$" - ^) 
2 4 A 



(39) 



at small and intermediate r. Polar contour plots of the lines of same energy density 
versus 9 are shown in Fig|2]and Fig|3]for the four new one-monopoles for values of 
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A = and A = 1. Hence from the polar contour plots of ED{r, 9) versus the angle 
6', we can conclude that these JEA type one-monopoles are axially symmetric 
about the z-axis. From these plots it is clearly seen that the magnitude of the 
Higgs magnetic field \r'^B^\ at infinity determines the shape of the one-monopole 
and hence a rough idea about the energy density distribution at small r. 

Near to the origin as r tends to zero, Bf — )■ and Bi ^ Bf . From the polar 
plots of \r^B^\ (Figjl| and the energy density contour plots (Fig|2]and Figjs]) of the 
four axially symmetric one-monopole, we are able to conclude that the (1, \,ki) 
one-monopole is stretched towards the positive ^-axis [lOj. The position of the 
monopole can be read from the zero of the Higgs field versus the x-z plane 3D plot 
(Figji]). When A = 0, the monopole is located along the z-ajds at z = 0.06 and 
when A = 1, the monopole is at z = 0.04. 

The one-monopole is however stretched in the opposite direction, that 

is, it is stretched along the negative 2;-axis. In this case the monopole is located 
along the z-axis at z = —0.25 when A = 0, and when A = 1 it is located at 
z = -0.32 (FigQ. 

However the (2, 1, ki) and (2, 3, ki) one-monopoles are of different shapes and 
both are equally compressed along both the positive and negative z-axis and are 
hence symmetrical about the z-axis and the x-y plane. Both these monopoles are 
located at the origin (Figj5|. 

The total energy ^ of these one-monopoles are also calculated numerically. 
Our numerical calculation gives the values as shown in Table [T] for 't Hooft- 
Polyakov one-monopole and the four axially symmetric one-monopoles for values 
of A = and A = 1. Since the minimum total energy cannot be less than one, 
we conclude that all the five one-monopoles have similar total energy equal to 
one. The deviation from 1 of the numerically calculated values gives the magni- 
tude of the numerical error involves in the calculation. As the complexity of the 
one-monopole gets higher, the numerical error gets bigger. 

When A = 1, the total energies of all the five monopoles are calculated to be 
1.29 in units of 47r^. Although all the monopoles seem to be having the same 
energies which is an invariant quantity for all gauge equivalent solutions, their 
energy densities do differ from one another as shown in Fig. 2 and 3. Since it is 
the energy density distribution in space that determines the shape of the monopole, 
we can conclude that the four new monopole are axially symmetric monopoles. 

We also notice that the position of the monopole is moved when the parameter 
p = I, Table |2} The monopole is located on the positive side of the 2;-axis when 
g = I and on the negative z-axis when q = ^■ 

5 Comments 

Although the total energy is equal to one for all the five monopoles when A = and 
equal to 1.29 when A = 1, the energy density distribution varies from one monopole 
to the next and the four new monopoles only possess cylindrical symmetries about 
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Table 1: The total energy in unit of An^ of the 't Hooft-Polyakov one-monopole 
and the four axially symmetric one-monopoles as calculated numerically when (a) 
A = and (b) when A = 1. 



{p, q, k) 


(0,0, A;) 




(1,1' ^i) 


(2,l,A;i) 


(2,3,A;i) 


(a) A = 


0.9996 


0.9994 


1.0138 


0.9982 


1.0122 


(b) A = 1 


1.2914 


1.2914 


1.2875 


1.2915 


1.2872 



Table 2: The position of the 't Hooft-Polyakov one-monopole and the four axially 
symmetric one-monopoles when (a) A = and (b) when A = 1. 



(P, k) 


(0,0, A;) 


(l,i^i) 


(1,1, A:i) 


(2,l,A;i) 


(2,3,A;i) 


(a) A = 


(0,0,0) 


(0,0,0.06) 


(0,0,-0.25) 


(0,0,0) 


(0,0,0) 


(b) A = 1 


(0,0,0) 


(0,0,0.04) 


(0,0,-0.32) 


(0,0,0) 


(0,0,0) 



the z-axis. We would also like to take note that our value for the energy when A = 1 
is quite close to the numerical value obtained by Bogomol'nyi and Marinov (1976), 
1.238, and the value obtained by Bais and Primack (1976), 1.30 [2]. However the 
numerical value obtained by Kleihaus et al. (2004) [5] is 1.41 which is quite 
different from our value of 1.29. 

We run the numerical solution of the 't Hooft-Polyakov one-monopole when 
(p, g, k) = (0, 0, k) with similar grid size of 70 x 60 as a pilot test for comparison 
with our results of the new one-monopoles. Hence we have shown that it is possible 
to construct a non radially symmetric one-monopole with finite energy. In this 
case, we have found axially symmetric one-monopole with finite energy. 

The gauge potentials of the Jacobi elliptic generalized solutions do not satisfied 
the radiation gauge condition. However even without this constraint, we can only 
find certain discrete configurations of axially symmetric one-monopole solutions. 
There do not exist a parameter that can be continuously varied from the 't Hooft- 
Polyakov monopole to the four various axially symmetric one-monopoles. The 
solutions exist for only certain discrete values of {p, q, k). The solution is not finite 
for other values of {p, q, k). At the moment we are also looking for numerical one- 
monopole solutions that can be continuously distorted from the 't Hooft-Polyakov 
one-monopole by varying a parameter in the the asymptotic large r exact solution. 
If God is willing, this work will be reported in the near future. 
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From Table [T| we also notice that for a fixed value of p, the complexity of the 
solutions and hence the numerical error increased with the next higher value of q. 

The uses of these one-monopole solutions are still not clear and why only certain 
shapes and sizes are allowed in Nature is also not yet understood. Hence further 
work can be done in this direction. 
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Energy Density, (p,q,k)=(1,Erf2,k1), A=0 Energy Density, (p,q,i()=(1,a2,i(1), J.=1 




Figure 2: A polar contour plot of the energy density of the {p,q,k) = (1, ^,ki) 
one-monopole and the {p,q,k) = (1, |, fci) one-monopole when (a) A = and (b) 
when A = 1. 
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(a) 



(h) 



Energy Density, (p,q,k)=[2,1,k1), J.=0 Energy Density, (p,q,k)=(2,1,k1), 1=1 




Energy Density, (p,q,k)=(2,3,k1), 1 = Energy Density, (p,q,k)=(2,3,k1), 1. = 1 




Figure 3: A polar contour plot of the energy density of the {p,q,k) = {2,1, ki) 
one-monopole and the {p,q,k) = (2,3, fci) one-monopole when (a) A = and (b) 
when A = 1. 
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Figure 4: A 3D plot of the Higgs field modulus |$| of the {p,q,k) = (1,^,^1) 
one-monopole and the {p,q,k) = (1, |, fci) one-monopole when (a) A = and (b) 
when A = 1. 
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Figure 5: A 3D plot of the Higgs field modulus |$| of the {p,q,k) = (2, 1,^1) 
one-monopole and the {p,q,k) = (2,3, fci) one-monopole when (a) A = and (b) 
when A = 1. 
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